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Abstract. Let R be a real closed field, Q C R[Yi, . . . ,Yi, Xi, . . . , X^], with 
degy(Q) < 2,degxiQ) < d,Q £ Q,#(Q) = m, and V C R[Xi,...,Xfe] 
with degjf(P) < d,P G 'P,#{'P) = s. Let S C R'^+'= be a semi-algebraic 
set defined by a Boolean formula without negations, with atoms P = 0, P > 
0, P < 0, P £ 7-" U Q. We describe an algorithm for computing the the 
Betti numbers of S generalizing a similar algorithm described in |6]. The 
complexity of the algorithm is bounded by {£smd)^'~"' ^ ' . The complexity of 
the algorithm interpolates between the doubly exponential time bounds for the 
known algorithms in the general case, and the polynomial complexity in case 
of semi-algebraic sets defined by few quadratic inequalities [6]. Moreover, for 
fixed m and k this algorithm has polynomial time complexity in the remaining 
parameters. 



1. Introduction and Main Results 

Let R be a real closed field and 5 C R'^ a semi-algebraic set defined by a Boolean 
formula with atoms of the form P > 0, P < 0, P = for P e C R[Xi, ...,Xk\. 
We call S a P-semi-algebraic set and the Boolean formula defining S a P-formula. 
If, instead, the Boolean formula has atoms of the form P = 0, P > 0, P < 0, P G V, 
and additionally contains no negation, then we will call S a P-closed semi-algebraic 
set, and the formula defining S a P-closed formula. Moreover, we call a P-closed 
semi-algebraic set S basic if the P-closed formula defining S* is a conjunction of 
atoms of the form P = 0, P > 0, P < 0, P eV. 

For any closed semi-algebraic set X C R*^, we denote by bi{X) the dimension of 
the Q- vector space, {X, Q) , which is the i-tli homology group of X with coeffi- 
cients in Q. We refer to [TOj for the definition of homology in the case of R being 
an arbitrary real closed field, not necessarily the field of real numbers. 

1.1. Brief History. Designing efficient algorithms of computing the Betti numbers 
of semi-algebraic sets is an important problem which has been considered by several 
authors. We give a brief description of the recent advances and direct the reader 
to the survey article [5] for a more detailed exposition. 

For general semi-algebraic sets the best known algorithm for computing all the 
Betti numbers is via triangulation using cylindrical algebraic decomposition (see 

1991 Mathematics Subject Classification. Primary 14P10, 14P25; Secondary 68W30. 

Key words and phrases. Betti numbers, Quadratic inequalities. Semi-algebraic sets. 

The first author was supported in part by an NSF grant CCF-0634907. The second author 
was supported by an NTU startup grant. Part of this work was done when the authors were 
visiting the Institute of Mathematics and Its Application, Minneapolis, IRMAR (Rennes) and the 
Nanyang Technological University, Singapore. 



2 



SAUGATA BASU, DMITRII V. PASECHNIK, AND MARIE-FRANgOISE ROY 



for example |10j ) whose complexity is doubly exponential in the mimber of vari- 
ables. There have been some small advances in obtaining singly exponential time 
algorithms for computing some of the Bctti numbers, but wc arc still very far from 
having an algorithm for computing all the Bctti numbers of a given scmi-algebraic 
set in singly exponential time. Singly exponential time algorithms for computing 
the number of connected components of a semi-algebraic set S (i.e. bo^S)) has 
been known for quite some time [HI [HI [T71 [2Ql [13] . More recently, an algorithm 
with singly exponential complexity is given in [11] for computing the first Betti 
number of semi-algebraic sets. The above result is generalized in |12j . where a 
singly exponential time algorithm is given for computing the first £ Betti numbers 
of semi-algebraic sets, where £ is allowed to be any constant. Finally, note that 
singly exponential time algorithm is also known for computing the Eulcr-Poincare 
characteristic (which is the alternating sum of Betti numbers) of semi-algebraic sets 

HI- 

In another direction, several researchers have considered a special class of semi- 
algebraic sets - namely, semi-algebraic sets defined using quadratic polynomials. 
While the topology of such sets can be arbitrarily complicated (since any semi- 
algebraic set can be defined as the image under a linear projection of a semi- 
algebraic defined by quadratic inequalities), it is possible to prove bounds on the 
Betti numbers of such sets which are polynomial in the number of variables and 
exponential in only the number of inequalities [3l [8] . (In contrast a semi-algebraic 
set defined by a single polynomial of degree > 2 can have exponentially large Betti 
numbers.) Polynomial time algorithms for testing emptiness of such sets (where the 
number of inequalities is fixed) were given in [21 119j . A polynomial time algorithm 
(without any restriction on the number of inequalities) is given in |ni(see also [7]) 
for computing a constant number of the top Bctti numbers of semi-algebraic sets 
defined by quadratic inequalities. If moreover the number of inequalities is fixed 
then the algorithm computes all the Betti numbers in polynomial time. More 
precisely, an algorithm is described which takes as input a semi-algebraic set, S, 
defined by Qi > 0, . . . , Qm > 0, where each Qi S R[Yi, . . . , Xg] has degree < 2, and 
computes the top p Betti numbers of S, bk-i{S), . . . , 6fc_j,(S'), in polynomial time. 
The complexity of the algorithm is X^f^o (™)^2°<"'"'''''"" ^ p^j. g^ed m, we obtain 
by letting p ~ £, an algorithm for computing all the Betti numbers of S whose 
complexity is ^•^ ' ' . 

The goal of this paper is to design an algorithm for computing the Betti numbers 
of semi-algebraic sets defined in terms of partly quadratic systems of polynomials 
whose complexity interpolates between the doubly exponential time bounds for the 
known algorithms in the general case, and the polynomial complexity in case of 
semi-algebraic sets defined by few quadratic inequalities. Our algorithm is partly 
based on techniques developed in [S] , where we prove a quantitative result bounding 
the Betti numbers of semi-algebraic sets defined by partly quadratic systems of 
polynomials. Before stating this result we introduce some notation that we are 
going to fix for the rest of the paper. 

Notation 1. We denote by 

• Q C R[Yi, . . . ,Yi, Xi, . . . , Xk], a family of polynomials with 
degy.(Q) < 2,deg;,(Q) < d,g G Q,#(Q) = m, 

and by 
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• V C R[v'Ci, . . . , Xk] a family of polynomials with 

degx{P)<d,P ^V,i^{V) = s. 

The following theorem is proved in [9]. 

Theorem 1.1. Let S C R^^*^ be a {V U Q)-closed semi- algebraic set. Then 

b{S) < £\Ois + e + m)£d)'^+2™. 

In particular, for m<£, we have b{S) < f{0{s + e)£d)''+^'" . 

The above theorem interpolates previously known bounds on the Betti numbers 
of general semi-algebraic sets (which are exponential in the number of variables) 
mi [U m [HI , and bounds on Betti numbers of semi-algebraic sets defined by 
quadratic inequalities (which are exponential only in the number of inequalities 
and polynomial in the number of variables) [3l . Indeed we recover these extreme 
cases by by setting £ and m (respectively, s, d and k) to 0(1) in the above bound. 

1.2. Main Results. The main result of this paper is algorithmic. We describe an 
algorithm (Algorithm[5]below) for computing all the Betti numbers of a closed semi- 
algebraic set defined by partly quadratic systems of polynomials. The complexity 
of this algorithm interpolates the complexity of the best known algorithms for 
computing the Betti numbers of general semi-algebraic sets on one hand, and those 
described by quadratic inequalities on the other. 

Definition 1.2 (Complexity). By complexity of an algorithm we will mean the 
number of arithmetic operations (including comparisons) performed by the algo- 
rithm in R. We refer the reader to [TUl Chapter 8] for a full discussion about the 
various measures of complexity. 

We prove the following theorem. 

Theorem 1.3. There exists an algorithm that takes as input the description of a 
{V U Q)-closed semi-algebraic set S (following the same notation as in Theorem 
and outputs its Betti numbers bo{S), . . . ,bi^i;-i{S) . The complexity of this 
algorithm is bounded by (£smd)'^ * ^ ' . 

The algorithm we describe is an adaptation of the algorithm in [5], to the case 
where there are parameters, and the degrees with respect to these parameters could 
be larger than two. In addition, in this paper we also treat the case of general VUQ- 
closed sets, not just basic closed ones as was done in [6]. We also provide more 
details and analyze the complexity of the algorithm more carefully, in order to take 
into account the dependence on the additional parameters. 

1.3. Significance from the point of view of computational complexity the- 
ory. The problem of computing the Betti numbers of semi-algebraic sets in general 
is a PSPACE-hard problem. We refer the reader to [6] and the references contained 
therein, for a detailed discussion of these hardness results. On the other hand, as 
shown in [5], the problem of computing the Betti numbers of semi-algebraic sets 
defined by a constant number of quadratic inequalities is solvable in polynomial 
time. This result depends critically on the quadratic dependence of the variables, 
as witnessed by the fact that the problem of computing the Betti numbers of a 
real algebraic variety defined by a single quartic equation is also PSPACE-hard. 
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We show in this paper that the problem of computing the Betti numbers of semi- 
algebraic sets defined by a constant number of polynomial inequalities is solvable in 
polynomial time, even if we allow a small (constant sized) subset of the variables to 
occur with degrees larger than two in the polynomials defining the given set. Note 
that such a result is not obtainable directly from the results in [6] by the naive 
method of replacing the monomials having degrees larger than two by a larger set 
of quadratic ones (introducing new variables and equations in the process). 

The rest of the paper is organized as follows. In Section [5] we describe some 
mathematical results concerning the topology of sets defined by quadratic inequal- 
ities. We often omit proofs if these appear elsewhere and just provide appropriate 
pointers to literature. In Section [3] we describe our algorithm for computing all the 
Bctti numbers of semi-algebraic sets defined by partly quadratic systems of poly- 
nomials and prove its correctness and complexity bounds, thus proving Theorem 

2. Topology of sets defined by partly quadratic systems of 

polynomials 

In this section we recall a construction described in [9j that will be important 
for the algorithm described later. We parametrize a construction introduced by 
Agrachev in jTj while studying the topology of sets defined by (purely) quadratic 
inequalities (that is without the parameters Xi, . . . ,Xk in our notation). However, 
we do not make any non-degeneracy assumptions on our polynomials, and we also 
avoid construction of Leray spectral sequences as done in [Tj . 

We first need to fix some notation. 

2.1. Mathematical Preliminaries. 

2.1.1. Some Notation. For all a G i? we define 

sign(a) = if a = 0, 
= 1 if a > 0, 
= -lifa<0. 

Let ^ be a finite subset of R[Xi, . . . , Xk]. A sign condition on A is an element of 
{0, 1,— l}-^. The realization of the sign condition cr, TZ(a,I&), is the basic semi- 
algebraic set 

{xeR' \ /\ sign(P(x)) = a(P)}. 

A weak sign condition on A is an element of {{0}, {0, 1}, {0, -1}}A The real- 
ization of the weak sign condition p, TZ{p, R'^), is the basic semi-algebraic set 

{x e a sign(^(^)) e PiP)}- 
PeA 

We often abbreviate 7?.((t, R*") by TZ{(j), and we denote by Sign(^) the set of 
realizable sign conditions Sign(yt) = {a E {0, 1, -1}-^ | 7^(cr) 7^ 0}. 

More generally, for any A C R[^i, . . . ,Xk] and a ^-formula <&, we denote by 
??.($, R'^), or simply 7^($), the semi-algebraic set defined by $ in R*^. 
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2.1.2. Use of Infinitesimals. Later in the paper, we extend the ground field R by 
infinitesimal elements. We denote by R(C) the real closed field of algebraic Puiseux 
series in C with coefficients in R (see HU] for more details). The sign of a Puiseux 
series in R(C) agrees with the sign of the coefficient of the lowest degree term in C. 
This induces a unique order on R(C) which makes ( infinitesimal: ^ is positive and 
smaller than any positive element of R. When a G R(C) is bounded from above 
and below by some elements of R, lim^(a) is the constant term of a, obtained by 
substituting for C, in a. We denote by R(Ci, • ■ • , Cn) the field R(Ci) • • • (Cn) and in 
this case Ci is positive and infinitesimally small compared to 1, and for 1 < z < n — 1, 
(^i+i is positive and infinitesimally small compared to Q, which we abbreviate by 
writing < Cn Ci 1- 

Let R' be a real closed field containing R. Given a semi-algebraic set S* in R'^, 
the extension of S to R', denoted Ext(5, R'), is the semi-algebraic subset of R'*^ 
defined by the same quantifier free formula that defines S. The set Ext(S', R') is 
well defined (i.e. it only depends on the set S and not on the quantifier free formula 
chosen to describe it). This is an easy consequence of the transfer principle (see for 
instance [10]). 

2.2. Homogeneous Case. 

Notation 2. We denote by 

• Q'* the family of polynomials obtained by homogenizing Q with respect to 
the variables Y, i.e. 

Q'' = {Q'' I QeQ}cR[Yo,...,Yi,Xi,...,Xkl 

where Q'' = Yo^Q{Yi/Yo, . . . , Ye/Yo, X^, . . . , X,,), 

• $ a formula defining a T'-closed semi-algebraic set V, 

• A'^ the semi-algebraic set 

(2.1) U {(f'^) I \y\ = ^ ^ Q(y,-^)<o A $(x)}, 

• the semi-algebraic set 

(2.2) W''= fl {{y,x) I I2/HI A g(y,x)<0 A $(x)}. 

Let 

(2.3) n = {uj eR"^ \\lu\ = l,Lu, <0,1 <i <m}. 

Let Q = {Qi,...,Qm} and Q'* = {Qi,...,Qt}- For w G wc denote by 
(w, Q'') € R[yo, ...,Ye,Xi,...,Xk] the polynomial defined by 

m 

(2.4) (c.,Q'')=^c..Qf. 

2 = 1 

For (w, x) G O X y, we denote by (oj, Q'^){-, x) the quadratic form in Yq, . . . ,Yg 
obtained from (w, Q^) by specializing Xi = Xi,l < i < k. 
Let B cnxS^ X V be the semi-algebraic set defined by 

(2.5) B^{{Lu,y,x)\u;en,yeS',xeV, {u, Q''){y,x) > 0}. 

We denote hy ipi : B ^ F and (p2 B ^ x V the two projection maps (see 
diagram below). 
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B 



F = nxv 





S'xV 



V 

The following key proposition was proved by Agrachev [1] in the unparametrized 
situation, but as we see below it works in the parametrized case as well. 

Proposition 2.1. The semi- algebraic set B is homotopy equivalent to A''. 

Proof. See [9]. □ 

We will use the following notation. 

Notation 3. For a quadratic form Q e R[ioj • ■ ■ j^^Jj we denote by index(Q) the 
number of negative eigenvalues of the symmetric matrix of the corresponding bi- 
linear form, i.e. of the matrix M such that Q{y) = {My,y) for all y e R^^^ (here 
(•,•) denotes the usual inner product). We also denote by Xi{Q),0 < i < i the 
eigenvalues of Q in non-decreasing order, i.e. 

Ao(0)<Ai(Q)<---<Af(Q). 

For F ~ X V a,s above we denote 

Fj^{{Lo,x)eF I inde^{{uj,Q''){;x))<j}. 

It is clear that each Fj is a closed semi-algebraic subset of F and we get a 
filtration of the space F given by 

Fo C Fi C • ■ • C = F. 

Lemma 2.2. The fibre of the map ipi over a point {uj,x) G Fj \ Fj^i has the 
homotopy type of a sphere of dimension £ — j ■ 

Proof See [g. □ 

For each G Fj \ Fj^i, let L'^{uj,x) C R*^^^ denote the sum of the non- 

negative eigenspaces of (cj, Q^){-, x). Since index((ti;, Q'')(-, x)) = j stays invariant 
as {uj,x) varies over Fj \ Fj-i, Lj'{uj,x) varies continuously with {uj,x). 

We denote by C the semi-algebraic set defined by the following. We first define 
for < j < ^ + 1 

(2.6) Cj^{iu,y,x) I iLu,x)eF,\Fj^i,yeL+{u;,x),\y\^l}, 
and finally we define 

e+i 

(2.7) C^IJC,. 

3=0 

The following proposition proved in [9] relates the homotopy type of B to that 
of C. 

Proposition 2.3. The semi- algebraic set C defined by (|2.7p is homotopy equivalent 
to B. 
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The following example which also appears in [9j illustrates Proposition 12.31 
Example 2.4. In this example m ^ 2, £ ^ 3, k = 0, and Q'' = {Q'l, Q^} with 

Q'i =Y^ + 2Y^ + Wl 

The set fi is the part of the unit circle in the third quadrant of the plane, 
and F = n \n this case (since fc = 0). In the following Figure [U we display 
the fibers of the map ^p^^iw) C B for a sequence of values of to starting from 
(—1,0) and ending at (0,-1). We also show the spheres, C fl ip^^{Ljj), of di- 
mensions 0,1, and 2, that these fibers retract to. At = (— 1, 0), it is easy to 
verify that index((a;, Q'*)) = 3, and the fiber (p^^(ijj) C i? is empty. Starting 
from u) = (— cos(arctan(l)), — sin(arctan(l))) we have index((a;, Q'')) = 2, and 
the fiber <p^^(a') consists of the union of two spherical caps, homotopy equiva- 
lent to S°. Starting from w = (— cos(arctan(l/2)), — sin(arctan(l/2))) we have 
index((a', Q'')) = 1, and the fiber ip'^^{w) is homotopy equivalent to S^. Finally, 
starting from = (— cos(arctan(l/3)), — sin(arctan(l/3))), index((ijj, Q'*)) = 0, 
and the fiber ip^^{uj) stays equal to to S^. 

^'OOOO 

Figure 1. Type change: ^ S° ^ ^ S^. is not shown. 
Let A e R[^i, . . . , ATi, . . . , Xfe, r] be the polynomial defined by 
A = det(T-Id,+i -M^.g.), 

= T^+i + c,r^ + • • • + Co, 

where Z ■ Q!" = YZi Z^Qt and each C, G R[Zi, . . . , Xi, . . . , Xu]. 

Note that for (oj, G x R*^, the polynomial A(a;, x, T), being the characteristic 
polynomial of a real symmetric matrix, has all its roots real. It then follows from 
Descartes' rule of signs (see for instance [IH]), that for each {lo,x) £ 17 x R'^, 
index((a;, Q'')(-, x)) is determined by the sign vector 

(s\gn{Ci{uj,x)), . . . ,sign(Co(w,a;))). 

More precisely, the number of sign variations in the sequence 

sign(Co(u;, x)), . . . , (-l)''sign(C,(w, x)), . . . , (-l)^sign(C^ (u;, x)), +1 

is equal to index((a;, Q'')(-, a;)). 
Hence, denoting 

(2.8) C — {Co, . . . , Cf} c R[Zi, . . . , Z„i, Xi, . . . , Xk], 

we have 

Lemma 2.5. Fj is the intersection of F with a C-closed semi- algebraic set Dj C 
R'"+'=, /or eac/i < j < £+ 1. □ 
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3. Computing the Betti numbers 

Wc now consider the algorithmic problem of computing all the Betti numbers of 
a semi-algebraic set defined by a partly quadratic system of polynomials. 

3.1. Summary of the main idea. The main idea behind the algorithm can be 
summarized as follows. 

By virtue of Proposition 12.11 in order to compute the Betti numbers of A'^ , it 
suffices to construct a cell complex, IC{B,V), whose associated space is homotopy 
equivalent to the set B defined by (|2.5[) . In order to do so, we first compute a 
semi-algebraic triangulation, h : A —>■ F, such that as {iu,x) varies over the image 
of any simplex cr G A, the index((u;, Q^){-, x)) stays fixed, and we have a continuous 
choice of an orthonormal basis, 

{eo(cr,w,x), . . . ,ee{(T,uj,x)} 

consisting of eigenvectors of the symmetric matrix associated to the quadratic form 

{tO,Q''){;x). 

Moreover, if index((a;, Q'^){-,x)) — j for {uj,x) G h{a), then (p-^ ^i^^,x) can be 
retracted to n span(ej {a,Lu,x), . . . , ee{(J, u!,x)), and the flag of subspaces deflned 
by the orthonormal basis, {eo(cr, w, cc), . . . , ei{a, uj, x)}, gives an efficient regular cell 
decomposition of the sphere S n span(ej((T, w, x), . . . , e^(tT, w, x)) into 2{£ — j + 1) 
cells, having two cells of each dimension from to £ — j (see Definition 13. 3p . 

Now consider a pair of simplices, cr, r G A, with cr ^ r. The orthonormal basis 
{eo(T, w, x), . . . , e^(T, w, x)}, deflned for {uj,x) G /i(r) might not have a continuous 
extension to h{a) on the boundary of /i(t). In particular, the cell decompositions 
of the fibers, S H span(ej(cr, x), . . . , e£(cr, x)), over points in (w, x) G h{a) 
might not be compatible with those over neighboring points in /i(r). In order to 
obtain a proper cell complex we need to compute a common refinement of the 
cell decomposition of the sphere over each point in {uj,x) G h{a) induced by the 
basis {eo(cr, w, x), . . . , e^(cr, w, x)}, and the one obtained as a limit of those over 
certain points in /i(t) converging to {uj,x). We need to further subdivide h{a) to 
ensure that over each cell of this sub-division the combinatorial type of the above 
refinements stays the same. Since a simplex cr G A can be incident on many 
other simplices of A, we might in the above procedure need to simultaneously 
refine cell decompositions of the sphere coming from many different simplices. In 
order to ensure (for complexity reasons) that we do not have to simultaneously 
refine cell decompositions coming from too many simplices, we thicken the simplices 
infinitesimally and as a result only need to refine at most m + k cell decompositions 
at a time. 

Before describing the construction of IC{B, V) in more detail, we need some 
preliminaries on triangulations. 

3.2. Triangulations. We first need to recall a fact from semi-algebraic geometry 
about triangulations of semi-algebraic sets, and then we define the notion of an 
Index Invariant Triangulation and give an algorithm for computing it. 

3.2.1. Triangulations of semi-algebraic sets. A triangulation of a closed and bound- 
ed semi-algebraic set S* is a simplicial complex A together with a semi-algebraic 
homeomorphism from | A| to 5. We always assume that the simplices in A are open. 
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Given such a triangulation we will often identify the simplices in A with their images 
in S under the given homeomorphism, and will refer to the triangulation by A. 

Given a triangulation A, the cohomology groups H*(S') are isomorphic to the 
simplicial cohomology groups H'(A) of the simplicial complex A and are in fact 
independent of the triangulation A (this fact is classical over R; see for instance 
[To] for a self-contained proof in the category of semi-algebraic sets). 

We call a triangulation hi : |Ai| — > 5 of a semi-algebraic set 5, to be a refinement 
of a triangulation /12 : IA2I — s- if for every simplex ai G Ai, there exists a simplex 
0-2 G A2 such that hi{ai) C /i2(o'2). 

Let 5*1 C ^2 be two compact semi-algebraic subsets of R'^ . We say that a semi- 
algebraic triangulation : | A| — + 6*2 of S'2, respects 5*1 if for every simplex cr G A, 
h{cr) n Si = h{(j) or 0. In this case, h~^{Si) is identified with a sub-complex of A 
and /i|h-i(Si) : h~^{Si) ^ is a semi-algebraic triangulation of Si. We will refer 
to this sub-complex by A | Si • 

We will need the following theorem which can be deduced from Section 9.2 in 
[11] (see also [10]). 

Theorem 3.1. Let Si C S2 C R'^ be closed and bounded semi-algebraic sets, and 
let hi : Ai Si,i — 1,2 be semi- algebraic triangulations of Si,S2- Then there 
exists a semi-algebraic triangulation h : A S2 of S2, such that A respects Si, A 
is a refinement of A2, and A\s-^ is a refinement 0/ Ai. 

Moreover, there exists an algorithm which computes such a triangulation with 
complexity bound (sd)'-^^^^ , where s is the number of polynomials used in the defi- 
nition of Si and S2, and d is a bound on their degrees. 

3.2.2. Parametrized eigenvector basis. Let M{u!,x) be the symmmetric matrix as- 
sociated to the quadratic form {uj, Q^){-,x) defined by (|2.4p . When M{ui,x) has 
simple eigenvalues for all possible choice of ui,x in some domain, there is a finite 
choice of orthonormal bases consisting of eigenvectors of M(a;,x). However, when 
Miio^x) has multiple eigenvalues, the number of choices of orthonormal basis of 
eigenvectors is infinite. In order to avoid the problem caused by the latter situation 
we are going to use an infinitesimal deformation as follows. 
Let < e <C 1 be an infinitesimal and 

(3.1) M^{uj,x) = (l-e)A/(u;,a;)-Hediag(0,l,2,...,£). 

Note that for every {uj,x) G x R*"' the eigenvalues of M^{io,x) in R(e) are dis- 
tinct and nonzero. Indeed, replace £ by < in the definition of Me(w,a:) and obtain 
Mt{uj,x). Observe that the statement is true if t = 1, since the matrix Mi{u!,x) 
has distinct eigenvalues. Thus, the set of t's in the algebraically closed field R[i] 
for which Mt{uj,x) has ^ -I- 1 distinct eigenvalues is non-empty, constructible and 
contains a open subset, since the condition of having distinct eigenvalues is a stable 
condition. Thus, there exists Eq > 0, such that for all t G (0, Eq), Mt{(-0, x) has € -I- 1 
distinct eigenvalues, and hence it is also the case for the infinitesimal e. 

Denote by A(Me(a;, x, T)) = det(r • Id^+i — M^{lo,x)) the characteristic poly- 
nomial of M^{iu, x). Let A C R[Zi, . . . , Zm, Xi, . . . , Xk] be a set of polynomials 
containing C (see (|2.8p ) and such that for every sign condition p e {0, 1, —1}'^ and 
every (a;,x) G TZ{p,D, x R*'), the Thom encodings of the roots of A(Me(a;, x), T) 
stay fixed, as well as the list of the non-singular minors of size £ in M^{uj, x,T) at 
each root of A{Me{iu, x), T). 
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Then choosing a non-vanishing minor and using Cramer's rule, we find {£ + 1)^ 
rational functions in the variables u, uj, x, T which give for every (u, w, x) G 'R-{p^ ^ x 
R*'"'""'^) the coordinates of an eigenvector v^{u,ijj,x,te) associated to the eigenvalue 

(where u denotes the co-ordinate left out in the non-singular ^ y, (. minor chosen 
for this eigenvalue in the application of Cramer's rule). We denote by ee(tJ,x,t£) 
the unit eigenvector cc, w, a;, when is an eigenvalue. 

If the eigenvalues A^^o < . . . < A^^i are in increasing order, we define 

ee.i{uj,x) = ee{uJ,X,\e^i). 
Note that for every {u},x) e 17 x R''' 

(lim(ee,o('^,a;)), . . . , lim(ee,f (w, x))) 

£ € 

is an orthonormal basis consisting of eigenvectors of M{uj,x). 

3.2.3. Index Invariant Triangulations. We now define a certain special kind of scmi- 
algebraic triangulation of F that will play an important role in our algorithm. 

Definition 3.2. (Index Invariant Triangulation) An index invariant triangulation 
of is a triangulation 

h: F 

of F, which respects all the realization of the weak sign conditions on V and A 
(see definition in I3.2.2[) . As a consequence, h respects the subsets Fj for ev- 
ery I C Q. Moreover, indcx((a', Q'')(-, x)), stays invariant as {u},x) varies over 
/i(cr), and the maps Ce^oicr) e^^i sending {lo^x) S h{<7) to the orthonormal 
basis eefl{w,x), . . . ,ee^i{uj,x), are uniformly defined. Note also that for every 
{uj,x) g h{<T), 

{ej (cr, uj,x)),..., ee{(T, u, x)) = {lim(eej (w, x)),..., \im{ee,e{uj, x))} 

is a basis for the linear subspace L^{uj,x) C R^^^, (which is the orthogonal com- 
plement to the sum of the eigenspaces corresponding to the first j eigenvalues of 

We now describe an algorithm for computing index invariant triangulations. 
Algorithm 1 (Index Invariant Triangulation). 

Input 

• A family of polynomials, Q'' = {g'i\ . . . , Qi^} C R[ro, . . . , Ff, Xi, . . . , X^], 
where each Q'' is homogeneous of degree 2 in the variables Yo, . . . ,Yi, and 
of degree at most d in Xi , . . . , Xk , 

• another family of polynomials, V C R[A'i, . . . , Xk], with deg(P) < d, P e 

• a P-closed formula $ defining a bounded P-closed semi-algebraic set V C 

R^ 

Output : an index invariant triangulation, 

h:A^F 

of F and for each simplex a of A, the rational functions eefl{a), . . . , ee,f((T). 
Procedure 
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Step 1. Let e > be an infinitesimal and let Z = (Zi, . . . , Zm)- Let be the 
symmetric matrix corresponding to the quadratic form (in Yq, . . . ,Y£) defined by 

Me{X, Z) = {1- e){ZiQ'l + ■■■ + ZrnQt) + eQ, 

where Q — X]i=o ■ Compute the polynomials 

(3.2) A(Z, X, T) = det(T • Id^+i - M^) = T^+i + C^T^ + • • • + Cq. 

Step 2. Using Algorithm 11.19 in [10] (Restricted Elimination), compute a family of 
polynomials A' C R[e] [Zi, . . . , Z^.Xi, . . . , X^] such that for each p e Sign(yl'), and 
{uj,x) G ^Z{p,^lxR'')nF the Thom encodings of the roots of A(w, a;,r) in R(e) and 
the number of non-negative roots of A(a;, x, T) in R(£) stay fixed, as well as the list 
of the non singular minors of size £ in M^{u!, x, T) at each root of A(Me(w, x), T). 
Let A C R[.^i, . . . , Zm,Xi, . . . , Xk] be the set of all coefficients of the polynomials 
in A' , when each of them is written as a polynomial in e. 

Step 3. Using the algorithm implicit in Theorem 13.11 (Triangulation) , compute a 
semi-algebraic triangulation, 

h:A^F, 

respecting all the realizations of the weak sign conditions on AUV. 
Step 4. For each simplex a of A, output the maps e£.o(o'): ■ • ■ i Se.e{o')- 

Complexity Analysis: The complexity of the algorithm is dominated by the 
complexity of Step 3, which is (s^md)^"'"^''' . □ 

Proof of Correctness: It follows from the fact that the triangulation respects 
all weak sign conditions on A that index((a;, Q'')(-, x)) is constant for (w, x) £ h{a) 
for any simplex ct of A, 

Since ee,o(o', w, x), . . ., ee,e{cr,uj,x) are orthonormal, so are eo(cr, w,x), . . ., 
ee{a,oj,x) for every {io,x) S h{a). Moreover, letting j = index((a', Q''}(-, x)) for 
{oJ,x) G /i(cr), we have that ee.j((T, x), . . . , e^^^ (cr, x) span the sum of the non- 
negative eigenspaces of Ms{lo,x), their images under the lim^ map will span the 
sum of the non-negative eigenspaces of M(a;, x). □ 

3.3. Computing Betti numbers in the homegeneous union case. Now that 
we obtained an Index Invariant Triangulation A, our next goal is to construct a cell 
complex JC{B,V) homotopy equivalent to B (see Notation l2.5|) that will be used to 
compute the Betti numbers of (see Notation . The cell complex IC{B,V) is 
obtained by glueing together certain regular cell complexes, JC{a), where cr G A. 

3.3.1. Definition o/C(A). Let 1 ^ Eq ^ ei ^ • • • ^ £m+k > be infinitesimals. 
For T e A, we denote by Dt the subset of f defined by 

Dr = {v E f I dist(u, 6) > edim(e) for all 9 -< ct}, 

where dist refers to the ordinary Euclidean distance. Now, let a ^ r be two 
simplices of A. We denote by Dc^^r the subset of f defined by 

Da,T = {v e f I dist(w, ct) < £dim(cr), and dist(u, 9) > £dim(9) for all 9 -< ct}. 

Note that 

|A| = IJ U y D^^r- 

o-eA cr,T6A,(T^r 
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Figure 2. The complex A. 




Figure 3. The correspondmg complex C(A). 



Also, observe that the various Dr's and D^r^s are all homeomorphic to closed 
balls, and moreover all non-empty intersections between them also have the same 
property. Thus, the union of the Dt^s and D^^s together with the non-empty 
intersections between them form a regular cell complex, C(A), whose underlying 
topological space is |A| (see Figures [5] and [3]). 

3.3.2. Definition of K.{a) and K,{a,T) where a,T are simplices of A. We now as- 
sociate to each (respectively, D^^r) a regular cell complex, /C(cr), (respectively, 
/C((T, r)) homotopy equivalent to ip^^{h{Da)) (respectively, Lp^'^{h[D„^r))- 
For each tr g A, and (a;,a;) £ h{a)^ the orthonormal basis 



{eo(cr,w,x)), . . . ,ei{(7,u),x)} 
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determines a complete flag of subspaces, T{<t,uj,x), consisting of 
F"(cr,w,a;) =0, 

F"'"((T, w, x) — span(ef (cr, tJ, x)), 

F^(cr, oj, x) =span(ef (cr, ti^, x), ef_i((T, x)), 

Definition 3.3. For < j < £, let c^((T,a;,a;) (respectively, c~(CT,a;,a;)) denote the 
(£ — j)-dimcnsional cell consisting of the intersection of the F^~^^^{ij,lo,x) with 
the unit hemisphere in R^"*"^ defined by 

{yes' I {y,e,{a,u,x))>Q} 
(respectively, {yeS^ | (y, (cr, x)) < 0} ). 

The regular cell complex lC{a) (as well as /C((t, r)) is defined as follows. 

For each v € |A| and cr € A, let v{(j) € \a\ denote the point of \a\ closest to v. 

The cells of /C(a) are 

I y e cf{(7,oj,x), {uj,x) e h{c)}, 

where index((w, Q'')(-, x)) < j < and c G C(A) is cither itself, or a cell 
contained in the boundary of 

Similarly, the cells of K.{<j,t) are 

{{y,^,x) I y e c^(cr,/i(w(cr))),-i; = hr'^{uj,x) e c}, 

where index((a;, Q'*)(-, x)) < j < ^, c G C(A) is cither Z^o-.t itself, or a cell contained 
in the boundary of D^ -r- 

3.3.3. Definition of ]C{D), where D is a cell o/C(A). Our next step is to obtain cel- 
lular subdivisions of each non-empty intersection amongst the spaces associated to 
the complexes constructed above, and thus obtain a regular cell complex, IC{B, V), 
whose associated space, \IC{B, V)\, will be shown to be homotopy equivalent to B 
(Proposition 13.61 below) . 

First notice that \IC{(j',t')\ (respectively, |A^(cr)|) has a non-empty intersection 
with |/C(cr, t)| only if D^/^t' (respectively, D„i) intersects D„t. 

Let D be some non-empty intersection amongst the Do-'s and ZJ^.t's, that is D 
is a cell of C(A). Then D C |t| for a unique simplex r G A, and 

D^D^,,rr\---r\D„^^rr\Dr, 

with ui < 02 < ■ ■ ■ < (Tp < CTp+i = T and p < m + k. 

For each z, 1 < i < p + 1, let {/o(o'i, v), . . . , fi{<Ji^ f )} denote a orthonormal basis 
of R^+^ where 



fj{(J,,v) = limej((T„/i(to(a,) + {I - t)v{<Ji))),Q < j < £, 
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and let T(ai, v) denote the corresponding flag, consisting of 

=0, 

F^(ai,v) =span(/£(cr^,i;)), 

We thus have p+1 different flags, 
and these give rise to p + 1 different regular cell decompositions of S*". 











T(ay.^) 
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Figure 4. The cell complex IC'{D,v). 

There is a unique smallest regular cell complex, IC'{D, v), that refines all these cell 
decompositions, whose cells are the following. Let L C R^'''^ be any j-dimensional 
linear subspace, < j < £ + 1, which is an intersection of linear subspaces 
Li, . . . , ip+i, where Li e J'io'i, v),l < i < p + I < m + k + 1. The elements of 
the flags, T(<j\^ w), . . . , T{(jpj^\^v) of dimensions j + 1, partition L into polyhedral 
cones of various dimensions. The intersections of these cones with S^, over all such 
subspaces L C R^"*"^, are the cells of K,'{D,v). Figure [J] illustrates the refinement 
described above in case of two flags in R'^. We denote by IC{D,v) the sub-complex 
of 1C'{D, v) consisting of only those cells included in i+(cri, h{v{ai))) n S^. 

We now triangulate h(D) using the algorithm implicit in Thcorem l3.1l (Triangu- 
lation), so that the combinatorial type of the arrangement of flags, 

T{ai,v), . . .,T{ap+i,v) 
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and hence the combinatorial type of the cell decomposition IC'{D, v), stays invariant 
over the image, hoiO), of each simplex, 0, of this triangulation. 

Note that in case the eigenvalues of M{h{v)) are all distinct, we have that 
J-{ai, v) = ■ ■ ■ = J-'{ap+i,v) since the vectors 

fo{cr.t,v),...,fi{cr.t,v) 

is an orthonromal basis of eigen- vectors with each fj(ai,v) uniquely defined upto 
sign. Thus the cell decompositions of induced by the flags J-{ai,v) ~ ■■ ■ — 
J-"((Tp+i,w) are identical to each other and hence also to ]C'{D,v). 

However, if the eigenvalues of AI{h{v)) are not all distinct then the refinement 
]C'{D,v) can be non-trivial. For example suppose we have that Xa{h{v)) = ■■ ■ = 
A^(/i(u)),0 < a < l3 < £. Then in general the sub-flags consisting of subspaces 
Fi+i_l3(ai,v) C •■• C F£+i_„(<Ti,u) and ^£+i_/3(<Tj, w) C ••• C F(+i^a{o-.j,v) will 
in general not coincide for i ^ j. 

In this case the combinatorial type of the reflnement 1C'{D,v) is determined by 
the dimensions of the intersections amongst the subspaces 

Ft+i^p{ai,v), . . . ,F£+i_a(cri,t;), I <i <p+l. 

The dimensions of these intersections are determined by the minimal linear depen- 
dencies amongst the vectors fi{aj,v), a < i < /3, 1 < j < p -I- 1, and these are in 
turn determined by the ranks of matrices with at most £ + 1 rows of the following 
form. The rows of the matrix consists of at most p + 1 blocks, with the j-th block 
of the shape fa{aj,v), . . . , fa^ (uj, v), where a < aj < j3. (Note that every row of 
the above matrix consists of rational functions evaluated at a single root \a{h{v)) 
of A(M(/i(i;)), r), and this root is common to all the rows. This fact is important 
since in order to perform algebraic computations on the entries of the matrix we 
need to eliminate just one variable corresponding to this single root.) 

Introducing an infinitesimal 5 such that 1;3>(5^£>0, we note that for each 
< J < 

fj{ai,h~^{uj,x)) = liTaee,j{h{Sv{ai) + (1 - S)v{ai))) 

O 

= lirn ej{at,h{tv{ai) + (1 - t)v{ai))). 

We consider all matrices with at most £ + 1 rows consisting of blocks of the 
shape, fa{(jj,v),...,faj{'Jj,v), with < a < < /3 < ^, < j < m -I- fc, and 
\a{h{v)) = •■•:= \p{h{v)). The number of such matrices is clearly bounded by 

^0(m+k) 

Using the uniform formula defining e^jlai) and Proposition 14.7 of [10| . and 
Algorithm 8.16 in [lOj (for computing determinants over an arbitrary domain), we 
compute a family of polynomials in R[.^i, . . . , Zm-, Xi, . . . , Xk] such that over each 
sign consition of this family the rank of the given matrix stays fixed. 

Let 

Ad C R[2'i, . . . , Zm, Xi, . . . , Xf.] 

be the union of all these sets of polynomials. 

The combinatorial type of the cell decomposition lC'{D,v) will stay invariant 
as {ll},x) varies over each connected component of any realizable sign condition on 
Ad C R[Zi , . . . , Zm , Xi , . . . , Xk] ■ 
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Given the degree bounds on the rational functions defining {ee,o(o'), ■ • ■ , ee,£(cr)}, 
{to, x) £ h(a), and the complexity of Algorithm 8.16 in |10| . it is clear that the num- 
ber and degrees of the polynomials in the family Ad arc bounded by {s£md)'^"^ ^ ' . 
We then use the algorithm implicit in Theorem 13.11 (Triangulation), with Ad as 
input, to obtain the required triangulation. 

The closures of the sets 

{iy,u;,x) I yecelCiD,h-'iu;,x)), iu;,x) e hihDie))} 

form a regular cell complex which we denote by 1C{D). 

The following proposition gives an upper bound on the size of the complex IC{D). 
We use the notation introduced in the previous paragraph. 

Proposition 3.4. For each {lu,x) G h{D), the number of cells in }C{D , {u , x)) 

is bounded by £0{m+k) ^ Moreover, the number of cells in the complex JC{D) is 
bounded by {s£md)'^ * ^ ' . 

Proof. The first part of the proposition follows from the fact that there arc at most 
{£ + ly^+^-'+i choices for the linear space L and the number of (j — 1) dimensional 
cells contained in L is bounded by 2™'*''"" (which is an upper bound on the number of 
full dimensional cells in an arrangement of at most m + k hyperplanes). The second 
part is a consequence of the complexity estimate in Theorem 13.11 (Triangulation) 
and the bounds on number and degrees of polynomials in the family Ad stated 
above. □ 

3.3.4. Definition of ]C{B,V). Note that there is a homeomorphism 

: \IC{a,,T)\n^^\h{D)) ^ \IC{D)\ 

which takes each cell of \IC{ai, r)| n ipi^ {h{D)) to a union of cells in K,{D). We use 
these homeomorphisms to glue the cell complexes lC{<7i, r) together to form the cell 
complex 1C{B,V). 

Definition 3.5. The complex JC{B,V) is the union of all the complexes K,{D) 
constructed above, where we use the maps iD,ai to make the obvious identifications. 
It is clear that /C(i?, V) so defined is a regular cell complex. 

We have 

Proposition 3.6. \K.{B^V)\ is homotopy equivalent to B. 

Proof. We have from Proposition l2.3l that the semi-algebraic set C C i? (sec ()2.7|) for 
definition) is homotopy equivalent to B. We now prove that \IC{B, V) \ is homotopy 
equivalent to C which will prove the proposition. 

Let X^+k = \I^{B, V)\ and for < j < m -I- fc - 1, let Xj = lim.^. Xj+i. 

It follows from an application of the Vietoris-Smale theorem [^Hj that for each 
j,0 < j < m + k — 1, Fixt{Xj, R(eo, . . . , £j)) is homotopy equivalent to ^j+i- Also, 
by construction of IC{B, V), wc have that Xq = linig^ \IC{B, V)\ — C, which proves 
the proposition. □ 

We also have 

Proposition 3.7. The number of cells in the cell complex 1C{B,V) is bounded by 
{s£md)^'^ '"^^^ . 

Proof. The proposition is a consequence of Proposition 13.41 and the fact that the 
number of cells in the complex C(A) is bounded by (s^toc?)^"'"^'"'' . □ 
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3.3.5. Algorithm for computing the Betti numbers in the homogeneous union case. 
We now describe formally an algorithm for computing the Betti numbers of A'^ 
using the complex IC(B, V) described above. 

Algorithm 2 (Betti numbers, homogeneous union case). 
Input 

• A family of polynomials Q'* C R[io, ■ ■ ■ ,Yt, ^i, ■ ■ ■ , ^k], homogeneous of 
degree 2 in the variables ¥„,..., Y«, degxiQ'') <d,Q^ Q!" , #(2'') = to, 

• another family of polynomials V C R[Ari, . . . , Xk], with dcgx{P) ^ d,P ^ 

r,*iv) = s, 

• a "P-closed formula <J?(a;) defining a bounded P-closed semi- algebraic set 

V c 

Output 

• a description of the cell complex IC{B, V), 

• the Betti numbers of A'^ where the semi-algebraic set A'^ is defined by 

A" = U {{y,x) I |y| = 1 A Qiy,x) < A <i>(x)}. 

Procedure 

Step 1. Call Algorithm [1] (Index Invariant Triangulation) with input ,V and <& 

and compute h and A. 
Step 2. Construct the cell complex C(A) (following its definition given in Section 

ESI). 

Step 3. For each cell D £ C(A), compute, using the algorithm implicit in Theorem 
13.11 (Triangulation). the cell complex IC{D). 

Step 4. Compute a description of IC{B, V), including the matrices corresponding to 
the differentials in the complex C,{IC{B,V)). 

Step 5. Compute the Betti numbers of the complex C,{IC{B,V)) using hnear alge- 
bra. 

Complexity Analysis: The complexity of the algorithm is (s^md)^"'"^'^' , using 
the complexity of Algorithm [TJ □ 

Proof of Correctness: The correctness of the algorithm is a consequence of 
the correctness of Algorithm [1] and Proposition 13.61 □ 

3.4. Computing Betti numbers in the homogeneous intersection case. 

3.4.1. Definition of IC{Bi, V). We now define a subcomplex of IC{B, V) correspond- 
ing to a subset / C [to] . 

We first extend a few definitions from Section [2] 

For each subset / C [to] , we denote by Q J' the subset of Q'' of polynomials with 
indices in / and by Clj the subset of 

n^{uje R™ I \uj\ = 1, < 0, 1 < i < to}, 

obtained by setting the coordinates corresponding to the elements of [to] \ / to 0. 
More precisely, 

f^/ = £ R™ I \uj\ = 1, < 0, for i e I, and for i e [to] \ /}. 
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Note that we have a natural inchision il/ '—>■ 0[„j] = f2. 

Similarly, we denote hy Fj C F = F^m], the set ftj x V, and denote by Bi C 
V the semi-algebraic set defined by 
Bi = {iu;,y,x) \ cj e ni,y e S',x eV, {cj , Q'') (y , x) > 0}. 

We denote by 1^1.7 : Bi Fj and ip2,i ■ Bj ^ x V the two projection maps. 
Now we define IC{Bi, V) for every / C [m]. 

Definition 3.8. The complex IC{Bi,V) is the union of all the complexes IC{D) 
in C{Aj) , where C(A/) is the subcomplcx of C(A) consisting of cells contained in 
Ai = h-\F,). 

Using proofs similar to the ones give for K,{B, V), we have 
Proposition 3.9. \K.{Bi,V)\ is homotopy equivalent to Bj. □ 
Algorithm 3 (Computing the collection of lC{Bi,V), Id [1 . . . , 771]). 

Input 

• Q!' = {Qj , . . . , Qt) C R[>^), . . . , , Xi, . . . , Xfc], where each is homo- 
geneous of degree 2 in the variables Y^^ . . . ^Yi, and of degree at most d in 
Xi,. . . , Xk, 

• PC R[Xi, . . . , Xfc], with deg(P) <d,P eV, 

• a P-closed formula <i>(x) defining a bounded P-closed semi- algebraic set 
V C 

Output 

• For each subset / C [m], a description of the cell complex IC{Bj, V). 

• For each / C J C [m], a homomorphism 

z/,j:C.(B/,y)^C.(Bj,l/) 

inducing the inclusion homomorphism : H*(_B/, V) H*(i3j, V). 
Procedure 

Step 1. CaU Algorithm [2] to compute /C(P, V). 

Step 2. Give a description of IC{Bj,V) for each 7 C [m] and compute the matrices 
corresponding to the differentials in the complex C,{IC{Bj, V)). 

Step 3. For I (Z J C [m] with compute the matrices for the homomorphisms of 
complexes, 

i/.j:C.(/c(P/,y))^c.(/c(Pj,y)) 

in the following way. 

The complex /C(P/, V) is a subcomplcx of IC{Bj, V) by construction. Compute 
the matrix for the inclusion homomorphism, 

tij : C.(/C(P/, V)) ^ C.(/C(Pj, V)) 

and output the matrix for the homomorphism. 

Complexity Analysis: The complexity of the algorithm is (s£md)^°''"^'"' , using 
the complexity of Algorithm [TJ □ 

Proof of Correctness: The correctness of the algorithm is a consequence of 
the correctness of Algorithm [1] and Proposition 13.61 □ 
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3.4.2. Algorithm for computing the Betti numbers in the homogeneous intersection 
case. Let M^'' C X R*-' be the semi-algebraic set defined by 

W''^ f] {{y, x) I |yH 1 A Qiy, x)<0 A ^x)}, 

using Notation [2] 
Then 

(3.3) R,{W'') = H*(Tot.(A/-.,.(/C(B,y)))), 
where N»^,{JC{B, V)) is the bi-complex 

(3.4) Np,,{lC{B,V))= C,(/C(i?j,y)), 

JC[m],#(,/)=p+l 

with the horizontal and vertical differentials defined as follows. The vertical differ- 
entials 

(3.5) dp,q : Np,,{lC{B,V)) ^ Np,q-i{K.{B,V)), 
are induced by the boundary homomorphisms, 

9g : C^{]C{Bi,V)) ^ C,_i(/C(i?,7,l/)), 
and the horizontal differentials 

(3.6) <5p,, : AAp,,(/C(B, V)) ^ A/'p+i,,(/C(B, V)) 
are defined by 

where J C [m], #(J) = _p + f , 

if&Up^,{K{B,V))= C.(/C(Bj,y)), 

JC[m],#(J)=p+l 

and for I C J C [m] 

: C.{IC{Bi,V)) ^ C,{IC{Bj,V)) 

denotes the homomorphism induced by inclusion. 
For a proof of (|3.3p see [5] . 

Using (|3.3[) , we are able to compute the Betti numbers of using only linear 
algebra, once we have computed the various complexes lC{Bi, V), as well as the ho- 
momorphisms ijj for all / C J C [to] using Algorithm[3l Moreover, the complexity 
of this algorithm is asymptotically the same as that of Algorithm [3l 

We now formally describe this algorithm. 

Algorithm 4 (Betti numbers, homogeneous intersection case). 
Input 

. A family of polynomials, Q'' = {Q?, . . . , Q^} C R[ro, . . . , F,, Xi, . . . , X^], 
homogeneous of degree 2 with respect to Yq, . . . ,Y(, degxiQ'^) < d,Q'^ £ 
Q^ 

• another family, V C R[Xi, . . . , AT^] with degj^ (P) < d,P e V ,4{V) ^ s, 

• a formula $ defining a bounded P-closed semi-algebraic set V. 
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Output the Betti numbers bi{W^), where W'^ is the semi-algebraic set defind by 
W''^ f] {(y,x) I |y| = 1 A QHv.^) < A 

Procedure 

Step 1. CaU Algorithm [3] (Computing the collection of IC{Bi,V)) to compute for 
each / C J C [m], the complex C,{IC{Bi ,V)) using the natural basis consisting of 
the cells of IC{Bi, V) of various dimensions, as well as the matrices in this basis for 
the inclusion homomorphisms 

^/.J : G.{IC{Bi, V)) ^ C.{IC{Bj, V)). 

Step 2. Using the data from the previous step, compute matrices corresponding to 
the differentials in the complex, Tot,(A/',,,(/C(B, F))), where Af,,,{IC{B,V)) is the 
bi-complex described by (|3.4|) - ()3.6p . 

Step 3. Compute, using linear algebra subroutines 

b,{W'') ^ H,((Tot.(A/-.,.(/C(B, V))))). 

Complexity Analysis: The complexity of the algorithm is dominated by the first 
step, whose complexity is (s£md)^°' ^ using the complexity of Algorithm[31 □ 

Proof of Correctness: The correctness of the algorithm is a consequence of 
the correctness of Algorithm [3] and (|3.3p . □ 

3.5. Computing Betti numbers of general V U Q-closed sets. Let S C R^'*''^ 
be a semi- algebraic set defined by a T' U Q closed formula <I>. 

Let Eq denote the set of all possible weak sign conditions on the family Q, i.e. 

Sq = {0,{0,1},{0,-1}}2. 

In the last section we defined a bi-complex J\f,^,{IC{B, V)) whose total complex 
has homology groups isomorphic to these of the semi-algebraic set = T^ip'^ H 0) , 
where p £ Eg is given by p{Qi) = {0,-1} for each i,l < i < m, and is 
obtained from p by replacing each Qi € Qhy Qi\ We now generalize this definition 
to the case of multiple weak sign conditions. More precisely, given a set E = 
{pi, . . . , ptv} C Eg, we define a corresponding bi-complex having properties similar 
to that of A/',,,(/C(i?, V)), but now with respect to E instead of a single weak sign 
condition p. 

Without loss of generality we can write <!> in the form 

Y pA<j)p, 

where each i^p is a P-closcd formula. 
Let 

Wp=7e(pA0p,R^+^), 
Let 1 ^ e > be an infinitesimal, and let 
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Po=e\xl + --- + Xl)-l, 
and Sb C R(£)^+'' be the semi- algebraic set defined by 

Sb = []{{v,x) I Qo(y) < A P^{x) < A $(2;)}. 

We denote by <i>f, (resp. 0p^f,) the formula {Qo{y) < 0) A {Po{x) < 0) A $ (resp. 
{Po{x) < 0) A 4>p-) 

Let S^, W^f, c X R(e}'' be the sets defined by $b and p A (Qg < 0) A (/)p,b 

respectively on x R(e)'" after replacing each S Q by in the formulas 
and p. 
Let 

Fp,, = 7e(0p,6,R(e)^). 
Let Q^^ = {±Q'' I e Q''}, and let 1C{B, V) denote the complex constructed by 
Algorithm [HI with input the families of polynomials Q^, Vt = V D {Pa}, and the 
semi-algebraic subset V = Bk{0, 1/e)- 

It follows from the correctness of Algorithm |4] that for each p g Sg, there exists 
Jp C Q5-, and a subcomplcx, IC{Bj^,Vp,b) C IC{B,V), such that 

H,(Tot.(AA...(/C(Bj^,yp,b))) - H,(VF;',,). 

More generally, for any S = {pi,...,p7v} C Sg, there exists a subcomplcx, 
IC{Bj^, Vp^^b<^ - ■ ■ VpN.b) C IC{B, V), such that the homology groups of the complex 

(3.7) Cs,. -Tot.(AA.,.(/C(Bj^,yp,.,n---yp„,,)) 

are naturally isomorphic to those of VFp'f,, where p is the common refinement of 
Pi, . . . , ptv defined by 

N 

(3.8) p(F) = f|p,(F) 

1=0 

for each P & A. Moreover, for E C S' C Sg, there exists a natural homomorphism, 

«s,s' : Cs'_, — *■ Cs,, 
such that the induced homomorphism, 

: H»(Cs',.) ^ H,(Cs,.) 
is the one induced by the inclusion 

n Kb n K^- 

Definition 3.10. Let C,($) denote the complex defined by 

(3.9) C.($) = Tot.(A/-.,.($)), 

where 

(3.10) KA^)^ Cs,,. 

ECX;a,#(E)=p+l 
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The vertical and horizontal homomorphisms in the complex A/',,,($) are induced 
by the the differentials in the individual complexes Cs,, and the inclusion homo- 
morphisms is,E' respectively. 

By the properties of the complexes Cs^, stated above and the exactness of the 
generalized Mayer- Vietoris sequence, we obtain 

Theorem 3.11. 

H,(5,^)-H4C.($)). 

We are now in a position to describe formally the algorithm for computing all 
the Betti numbers of a given V U Q-closed set S. 

3.5.1. Description of the algorithm in the general case. 

Algorithm 5 (Betti numbers, general case). 
Input 

• A family of polynomials Q = {Qi, . . . , Q„^} C R[Fi, . . . ,Yi, Xi, . . . , Xk], 
with AegyiQi) < 2, deg^(Q,;) <d,l<i<£, 

• another family of polynomials V C R[^i, . . . , Xk] with deg(P) < d,P & P, 

• a Q U "P-closed semi- algebraic set S defined by a Q U 'P-closed formula 
Output the Betti numbers feo(<S), . . . , bk+e-i{S). 

Procedure 

Step 1. Define Qo = e'^iY^^ + . . . + F/) - 1, Pq = eliXf + ... + Xl)-l. Replace 5* 

by 7^(5',R(e)) n (7^(Po < 0) x n{Qo < 0)). 
Step 2. Define Q'l = {±Q'' \ G Q''} U {Qg}, and let ]C{B,V) denote the 

complex constructed by Algorithm^ with input the families of polynomials Q\,Vb, 

and the semi-algebraic set V = Pfc(0, C R(£)*'. 
Step 3. Compute, using the definitions given above, the matrices corresponding to 

the differentials in the complex C,($). 
Step 4. Compute, using linear algebra subroutines, for each i,0<i<fc + £ — 1 

hiS'i) = H,(C.(<J>)). 

Step 5. Output for each i,0<i<fc + ^— 1, 

h{s) = h,{s'i). 

Proof of Correctness: The correctness of the algorithm follows from Theorem 
13.111 and the correctness of Algorithm [31 □ 

Complexity Analysis: Since #(Sg) = 3™, the number of subsets that enters 
in the definition of A/',.,($) (cf. p.lOp ') is at most 2^ . The complexity of the 
algorithm is now seen to be (simd)'^'^'' ^ ' , using the complexity of Algorithm [31 □ 

Proof of Theorem ] The correctness and complexity analysis of Agorithm[5lalso 
proves Theorem ll.3l □ 
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